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ITERATIVE DECODING AND EQUALIZING METHOD FOR HIGH SPEED 
COMMUNICATIONS ON MULTIPLE ANTENNA CHANNELS DURING 
TRANSMISSION AND RECEPTION 

The present invention relates to the field of 
digital communication . 

It relates more particularly to an iterative 
decoding and equalizing method and device for high-speed 
communication over multiple transmit antenna channels and 
more particularly over frequency-selective channels with 
multiple transmit and receive antennas. 

GENERAL TECHNICAL FIELD 

With the development of antenna technologies, 
communications models based on TDMA, CDMA, POFDM and 
combinations thereof are being systematically revised and 
expanded to encompass MIMO channels. Two major space- 
time coding classes are envisaged, which have different 
philosophies : 

BLAST transmission developed by Bell Labs aims to 
use space multiplexing to increase the bit rate that can 
be transmitted over the channel. There are two ways to 
combine BLAST with error corrector channel coding, 
yielding two different layered space-time coding models. 
In the first (ST BICM) model, a single coding operation 
linking the various layers is applied to the data before 
space- time interleaving and space demultiplexing [1, LST- 
I, p. 1350] . In the second (MAC) model, the data is 
space demultiplexed and then coded and interleaved 
independently on each layer [1, LST-II, p. 1350] . 

Diversity schemes use space diversity to combat the 
effects of channel fading more effectively and to improve 
performance. There are STTC schemes, which improve 
coding, and STBC schemes (and their generalizations OD, 
LDC, etc.), which do not lead to coding gain 5 and are 
often combined with the best trellis coded modulations 
known in the art for the Gaussian channel. 



Recent theoretical contributions have proved that 
improvements in space multiplexing and space diversity- 
are linked by a relationship of compromise. For example, 
the two above-mentioned approaches cannot increase 
5 capacity in the same proportions for arbitrarily high 
numbers M, N of transmit antennas and receive antennas . 
It can be shown that transmit diversity schemes are 
optimized from the information theory point of view only 
in the MISO configuration {M, N) = {2,1}, which seriously 
10 calls into question their relevance to high bit rate 

communication over MIMO channels using a large number of 
antennas . 

Only the ST-BICM approach guarantees linear 
augmentation of capacity with min {M, N} . 

15 This is why it is necessary to opt for this type of 

transmission and to concentrate design effort on a 
suitable receiver guaranteeing performance as close as 
possible to the fundamental limits. Unlike STBC, ST-BICM 
does not ensure mutual orthogonality of the data streams 

20 transmitted by the various antennas. 

In the case of transmission over a frequency- 
selective MIMO channel, this coding strategy therefore 
calls for conjoint detection of data at the receiver to 
eliminate multiple access interference (MAI) in the space 

25 domain and intersymbol interference (ISI) in the time 
domain . 

The conjoint data detection function constitutes 
the delicate and complex element of the receiver, 
especially as ST-BICM inherently necessitates iterative 
30 processing between channel detection and decoding to 
capture all of the diversity in the space domain. 



GENERAL DESCRIPTION OF THE PRIOR ART 

Iterative decoding on non-frequency- selective MIMO 

channels 

Iterative decoding of ST-BICM on MIMO channels that 
5 are not frequency selective constitute a special case of 
the problem to be solved under the two-fold assumption of 
perfect ISI cancellation and recovery of all the energy 
available for each symbol. 

The problem may then be restated in the following 
10 terms: ''cancellation of MAI on a flat MIMO channel in the 
presence of spatially colored noise". It necessitates 
calculation at all times of the APP for the bits of the 
symbols of the vectors transmitted (statistics on the 
data to be forwarded to the decoder) . 
15 Two methods are described in the literature: 

- Exact calculation of the APP based on an 
exhaustive list (MAP) . The complexity of this method is 
prohibitive, since it increases exponentially with the 
number of transmit antennas M and is a polynomial 

20 function of the order Q of the constellations employed. 

- Approximate calculation of the APP from a reduced 
list of candidate vectors generated by a sphere decoding 
algorithm [2] . The sphere decoding algorithm has at 
least two advantages: it preserves the MAP criterion and 

25 has a complexity in O (M3) that is insensitive to the 
order Q of the constellations employed. Note that the 
sphere decoding algorithm may be regarded as a variant of 
the Fano sequential decoding algorithm. 

- Approximate calculation of the APP from a short 

3 0 list of candidate vectors generated by a sequential stack 
algorithm. 

The sphere decoding algorithm has been intensively 
studied for MIMO channels that are not frequency- 
selective. Generalizing to the situation of frequency- 
35 selective MIMO channels is non-trivial. 
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A brutal approach would be block replication. This 
approach is of little interest as it artificially 
increases the size of the search space (and therefore 
increases complexity) and introduces redundant decoding 
5 of symbols • 

The algorithm described in [2] is derived for MIMO 
channels that are not frequency selective. The jury is 
still out on generalizing it to the situation of 
frequency-selective MIMO channels, 

10 The performance of these two algorithms is 

essentially determined by the quality of the list of 
candidates generated, which is preferably of fixed size. 
It is crucial for it to include the best candidate 
decoding a 1 bit and the best candidate decoding a 0 bit 

15 in any position. 

Iterative decoding on frequency-selective MIMO channels 

Several types of detectors inserted into iterative 
structures have been proposed. 

20 - Detectors with weighted inputs and outputs 

applying the MAP criterion based on the BCJR algorithm. 
The complexity of this receiver in O (QLxM) is clearly 
prohibitive for high order modulation and MIMO channels 
with a large number of inputs and high memory demand. 

25 - Sub-optimal weighted input and output detectors 

based on efficient search algorithms in trellises greatly 
reduced in terms of number of states [3] and preceded, 
where appropriate, by minimum phase filters. This 
approach is still limited by its complexity in O (QLrxM) , 

3 0 where Lr represents the reduced constraint length. In 

particular, the use of constellations with a large number 
of states is excluded. 

Transposing non-linear detectors constructed from 
iterative linear cancellation by interference (see [4] in 

35 a multi-user detection context) to the MIMO context is 
generally based on the analogy between different users 
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and antennas . Their implementation therefore requires 
one filter per antenna and their performance is limited 
by the non- contiguous MMSE per antenna approach, which 
has the advantage that its complexity is a polynomial 
5 function of the number of send antennas, 

SUMMARY OF THE INVENTION 

An object of the invention is to propose an 
advanced receiver for high bit rate digital transmission 

10 over f rec[uency-selective channels with multiple transmit 
cind receive antennas that is not very complex - and in 
particular does not necessitate large computation 
powers - whilst providing for processing interference 
between antennas in the space domain and interference 

15 between symbols in the time domain. 

To this end, the invention proposes an iterative 
decoding and equalizing device for high bit rate 
communication over frequency-selective channels with 
multiple transmit and receive antennas, said device 

20 comprising a decision feedback equalizer adapted to 
receive data from different receive antennas and 
including a forward filter and a recursive backward 
filter fed with calculated weighted reconstituted data 
from the output of a decoder fed by decision means and 

25 the device further including means for subtracting the 
output of said backward filter from the output data of 
the forward filter, whereby the subtracted data is fed to 
the input of the decision means with the output of the 
decoder and the decision means produce a statistic which 

3 0 is forwarded to a channel decoder with weighted inputs 

and outputs and said decision means take into account the 
space noise correlation at the output of the subtraction 
means and the decision means and the decoder are 
separated by space-time interleaving at a binary level, 

3 5 which device is characterized in that the forward filter 
and the backward filter are iteratively adapted to 



minimize the mean square error at the output of the 
subtracter. 

A device of the above kind advantageously has the 
following complementary features, individually or in all 
5 technically possible combinations: 

- the decision means at the output of the 
subtraction means of the equalizer are of the space 
whitening type and followed by a sphere decoder; 

- the decision means at the output of the 

10 subtraction means of the equalizer are of the serial 

and/or parallel type (SIC/PIC) adapted to cancel residual 
space interference at the output of the subtraction means 
of the equalizer; 

- the space whitening is effected at the output of 
15 the subtraction means of the equalizer; 

- the space whitening is effected by the decision 

means ; 

- the space whitening is effected by the forward 
filter and the backward filter; 

20 - starting from a certain iteration, the forward 

filter is an adapted filter, the backward filter the Scime 
adapted filter minus the central coefficient. 

The invention stems from a novel approach that is 
radically different from those that consider the antennas 

25 as different users. The signal transmitted is seen as a 
T-dimensional modulation, where T is the number of 
transmit antennas. This changed point of view has 
considerable impact on the design of the receiver, which 
consists of a linear iterative equalizer considering at 

3 0 its input T-dimensional vector modulation convoluted with 
a frequency-selective channel, followed by a 
T-dimensional modulation detector (for example a sphere 
decoder) capable of generating flexible information for 
the channel decoder. This approach has two advantages 

35 over the prior art in that it is less complex, as it 
necessitates a single (vectorial) filter, and offers 
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better performance because it allows a choice of the 
vector detection criterion. 

The invention also relates to a system for high bit 
rate communication over frequency-selective channels with 
5 multiple transmit and receive antennas characterized in 
that it includes a receiver that includes an equalization 
and decoding device according of the above type. 

The transmit means are of the ST-BICM type, for 
example, which is advantageous. 

10 The invention further proposes an iterative 

decoding and equalizing method for high bit rate 
communication over frequency- selective channels with 
multiple transmit and receive antennas, using a decision 
feedback equalizer adapted to receive data from different 

15 receiving antennas and including a forward filter and a 
recursive backward filter fed with calculated weighted 
reconstituted data from the output of a decoder fed by 
decision means and using means for subtracting the output 
of said backward filter from the output data of the 

20 forward filter whereby the subtracted data is fed to the 
input of the decision means with the output of the 
decoder and the decision means produce a statistic which 
is forwarded to a channel decoder with weighted inputs 
and outputs, and said decision means take into account 

25 the space noise correlation at the output of the 

subtraction means and the decision means and the decoder 
are separated by space-time interleaving at a binary- 
level, which method is characterized in that the forward 
filter and the backward filter are iteratively adapted to 

3 0 minimize the mean square error at the output of the 
subtracter. 

DESCRIPTION OF THE DRAWINGS 

Other features and advantages of the invention 
3 5 emerge from the following description, which is purely 
illustrative and is not limiting on the invention, and 



should be read in conjunction with the appended drawings, 
in which: 

- Figure 1 shows the VBLAST concept; 

- Figure 2 shows a general model of ST-BICM 
5 communication; 

- Figure 3 shows the architecture of an iterative 
receiver conforming to one embodiment of the invention; 
and 

- Figure 4 shows the architecture of an iterative 
10 receiver that can be used from the second iteration, for 

example, or from subsequent iterations. 

DESCRIPTION OF ONE OR MORE EMBODIMENTS OF THE INVENTION 
BLAST technique and ST-BICM coding 

15 Figure 1 is a general diagram of a BLAST 

architecture . 

After space-time interleaving and space 
demultiplexing, data coded in a single vectorial encoding 
device EV is transmitted by a plurality of transmit 

20 antennas TX; at the other end of the MIMO channels, the 
transmitted signals are received by a plurality of 
receiver antennas RX that forward them to decoding means 
DE (MIMO channels) at the output whereof the data is 
recovered . 

25 The following description relates to frequency- 

selective MIMO channels. The presence of inter-symbol 
interference increases the complexity of the receiver. 

Figure 2 shows the ST-BICM general communications 

model . 

3 0 On transmission, the data is corrector channel 

coded (1) (convolutional code, turbocode, LDPC code, 
etc.), interleaved (2) at the binary level, space 
demultiplexed (3) and, for each layer, modulated (4-1 to 
4-M) . After shaping filtering (5-1 to 5-M) , the 

35 modulated data is forwarded to transmit antennas 6-1 to 
6-M. 
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Description of transmit and receive processing 

On transmission, the data is subjected to ST-BICM 

processing. 

This involves the following steps: 
5 - receiving a digital data stream at a given bit 

rate; 

- applying channel corrector coding 1 to generate a 
coded digital data stream; 

- interleaving (2) the coded digital data by means 
10 of a bit level interleaver; 

- space demultiplexing (3) the interleaved coded 
digital data stream to create a plurality of separate 
coded digital data streams called transmitter layers (M 
distinct streams or transmitter layers) . 

15 Then, for each of the M distinct coded digital data 

streams : 

- modulating (4-1 to 4-M) the interleaved coded 
digital data stream 15 in accordance with a modulation 
scheme to obtain a stream of modulated symbols; 

20 - filtering (5-1 to 5-M) the modulated stream; 

- transmitting (6-1 to 6-M) the modulated stream 
via its own transmitter antenna. 

Figure 3 shows the receiver of the transmission 
system. 

25 It includes a decision feedback equalizer defined 

by a forward filter 9, a subtracter 10, a decision 
algorithm 11 and a recursive backward filter 12. 

The recursive backward filter 12 is fed with 
weighted reconstituted calculated data from the output of 

30 a decoder 13 fed by the decision algorithm 11. 

The forward filter 9 and the backward filter 12 are 
determined iteratively to minimize the mean square error 
(MSE) of the MIMO equalizer, i.e. to minimize the error 
at the output of the subtracter 10. 

35 To this end they use an initial estimate 8 of the 

MIMO channel. 
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The processing employed is of the type described in 
appendix II, for example. 

Note that, in the context of the above processing, 
a vectorial estimate of the residual error and noise is 
5 used, the forward and backward filters being calculated 
block by block to minimize this vector. 

This vectorial processing simplifies calculation. 

The input of the decision means 11 receives data 
from the output of the subtracter 10 and from the output 
10 of the decoder 13. 

The algorithm may be of a different type, in 
particular of the SIC/PIC type (serial and/or parallel 
cancellation of residual space interference at the output 
of the subtracter means of the equalizer - see 
15 appendix III) . 

Alternatively, the algorithm may be a sphere 
decoder algorithm. 

An algorithm of this kind has a complexity in O (M^ ) 
(where M is the number of transmit antennas that is 
20 independent of Q. 

This enables the use of modulation with a large 
number of states with a view to increasing the bit rate. 

The output of the decision means 11 is forwarded to 
a space-time de-interleaving process 14 implemented at 
25 the binary level between said decision means 11 and the 
decoder 13 . 

The output of the decoder 13 is a bit probability. 
This probability is forwarded firstly to the input 
of the decision algorithm 11 and secondly, after space- 
3 0 time interleaving (15) , for weighted data reconstitution 
processing 16. 

The weighted data reconstituted in this way is 
forwarded to the input of the backward filter 12. 

The output of the backward filter 12 constitutes 
3 5 weighted regenerated data. 
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Furthermore, the error corresponding to the 
residual interference and noise that is injected at the 
input of the subtractor 10 is colored both in the time 
domain and in the space domain. 
5 Although the time correlation has little impact on 

the processing carried out by the equalizer, the space 
correlation has an essential role. 

This is why, in one embodiment, space whitening is 
effected by means of Cholesky factorization. 
10 This space whitening is advantageously effected at 

the output of the subtractor 10. 

In one particular embodiment, it may be effected by 
the sphere decoding algorithm 11 itself. 

It may equally be integrated into the forward and 
15 backward filters 9 and 12. 

Finally, in one embodiment, from a certain given 
number of iterations, for example from the second 
iteration, the forward filter 9 is advantageously 
replaced by an adapted filter. 
20 As the backward filter is deduced directly from the 

forward filter, convolution of the adapted filter with 
the channel minus the central coefficient therefore 
results in this case. 

This is shown in Figure 4, in which the forward 
25 filter 9 is an adapted filter MF while the backward 

filter 12 Bmf is the adapted filter convoluted with the 
channel minus the central coefficient . 
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APPENDIX I 

LIST OF ABBREVIATIONS 





APP : 


A Posteriori Probability 


5 


BCJR : 


Bahl, Cocke, Jelinek, Raviv (algorithm) 




BLAST : 


Bell Labs Layered Space Time 




CDMA : 


Code Division Multiple Access 




IC : 


Interference Cancellation 




ISI : 


Inter Symbol Interference 


10 


LDPC 


Low Density Parity Check 




MAI 


Multiple Access Interference 




MAP 


. Maximum A Posteriori 




MF 


; Matched Filter 




MIMO 


: Multiple Input Multiple Output 


15 


MMSE 


: Minimum Mean Square Error 




OFDM 


: Orthogonal Frequency Division Multiplex 




PIC 


: Parallel Interference Cancellation 




SIC 


: Serial Interference Cancellation 




SISO 


: Single Input Single Output 


20 


Soft-IC 


: Soft Interference Cancellation 




STBC 


: Space-Time Block Codes 




ST-BICM 


: Space-Time Bit -Interleaved Coded Modulation 




STTC 


: Space-Time Trellis Codes 




TDMA 


: Time Division Multiple Access 


25 


WMF 


: Whitening Matched Filter 
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APPENDIX II 
COMMUNICATIONS MODEL 

Consider a P-block channel with multiple inputs and 
5 outputs, T frequency-selective transmit antennas and R 
frequency-selective receive antennas of memory M, 

A. Space- time binairy interleaving coded modulation 

Let C denote a linear code of length Nc and of yield 

10 pc on ¥2 admitting a data vector u e F2^^'' and producing a 
coded word c e F2^. It is assumed that the yield pc 
includes tail bits if conventional codes are employed. 
The coded word enters a well-designed binary interleaver 
n, the output matrix A e F2'^'^*^^ whereof is segmented into 

15 P matrices A^ e F2'^'^*^^ , p = 1, P. The columns of the 
matrix A^ are vectors a^ [n] g F2'^'^, n = 1, L, called 
''symbol label vectors", containing T sub-vectors a^t [n] g 
F2*^, t = 1, T (one per input channel), with the stacked 
binary components aP<t,j>[n], ... aP<t,j>[n], where <t,j> 

20 represents the index (t - 1) q+ j . In each matrix A^, all 
the vectors a^[n] are modulated by a memoryless 
D-dimensional modulator over a signal set A c c° of 
cardinality |a| = 2*^ via bit coding with symbol m {0, 
1}"^ — A. Without loss of generality, it is assumed that 

25 A is a simple product of T complex constellations C of 
the same order 2^ and that Gray coding is used for each 
constellation. After transformation of the signal, this 
coded modulation process may also be considered as a 
space- time coded modulation scheme in which all the code 

3 0 word complex symbols PxLxt are divided into P matrices 

Xp e C^''^, p = 1, P, whereof the columns x^ [n] g c"^, 

n = 1, L are called ''constellation symbol vectors". 

The ^'symbol label vector" a^ [n] may be obtained from x^ 
[n] by simple inversion of bit coding with symbol 

35 Falling into the general class of space-time codes, this 
architecture is called space-time bit-interleaved coded 
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modulation (STBICM) [3] . From the Nyquist limited band 
ideal filtering assumption, the transmission yield (in 
bits per channel used) is: 

5 Figure 1 shows a diagram of this communications 

model . 

B. MIMO channel 

Let e c^^Tx(M+i) ^j^^ MIMO channel for the block 

10 p, where p = 1, . . P and H = {H^} the set of all the 

channels. Also let g c'^''^ and g C^""^ respectively be 
the "constellation symbol matrix" and the ''channel output 
matrix" . The vectorial channel output in the equivalent 
baseband at discrete time y^ [n] e at time n = 1, 

15 L may be written as follows: 

[n\ ^ {H {z)\ ^ [n] + [«) = ^ IF [k] x? (n - ib] 4- [n] (2) 

In the above equation, x^ [n] g C*^ are the 
''constellation symbol vectors" transmitted at the times 
n, the energy of each component being equal to Cx^/ [k] 

20 e C^'^'^ is the matrix of the number of taps k of the 
impulse response of the channel, w^ [n] e is the 
additive complex noise vector. The additive complex 
noise vectors w^ [n] are assumed to have a null, 
independent, identically distributed mean of complex 

25 Gaussian type with circular symmetry, and therefore 
follow the pdf N (0, o^I) . The channel H^, constant 
throughout the corresponding block, has a finite impulse 
response (FIR) of length M+1, the space symbol taps 
[0], . . . , [M] whereof are the random complex matrices 

3 0 RxT, of null mean, and the average power whereof 
satisfies the normalization constraints: 

E l^diag jx^H*" (fc] Iff (frj^l j - Tl (3) 
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in the case of an equal power system. The operator '*t" 
corresponds to the conjugate transposition operator. 
Equal channel memory for all the RxT possible links is a 
reasonable assumption given that the number of individual 
5 multichannel components is dictated predominantly by 
large structures and reflecting objects. 

MIMO interference cancellation iterative block 

10 A. Theory and scheduling 

The MIMO interference canceller iterative block 
processes each received data block Y^, p = 1, P 
separately using random probabilistic information fed 
back by the output decoder. All the signals and elements 

15 described are therefore indexed block by block. This 
dependency is sometimes dispensed with to simplify the 
notation. During each iteration I, a linear forward 
filter F"^ applied to each received symbol of vector y[n] 
produces the vector signal y^ [n] . Then an appropriately 

20 constructed estimate e^[n] of the MAI and of the ISI 

degrading x[n] is subtracted from the vector signal y''*[n] 
to produce [n] : 

z*H = y'H-e'tn] (8) 

The estimate of the vector e-'-En] comes from the 
25 output of the backward filter B-"- excited by the flexible 
decision attempt vector x^[n] on the vector x[n] , given 
the prior knowledge available (extrinsic probability 
distribution) on the iteration 1 - 1 . 

The description now turns to the core of the MIMO 
30 turbo -equalizer, namely the derivation of the MIMO 

forward and backward finite impulse response filters. 
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B. Calculation of forward and backward filters 

The description begins with a few basic 
manipulations of the instantaneous output y^[n] of the 
forward filter F^: 



Since: 



Af 



(22) 



(23) 



each sample y^[n] may be expanded as follows: 

y'H= E J2^mm^[n-i-k]+ ^ F-IilwIn-iJ (24) 

10 An equivalent matrix notation for the convolution 

of two filters He and is: 

[«] = F'H^ [n] + Ww^ [n] (25) 

in which: 

F' = [F' . ■ F^ [0] . . . F^ [Li^\] € C^x«^' (26) 

15 is the forward filter of order Lp = Lpi + Lf2 +1, where: 
H[0] H[l] ^-^ H[M] 

H[0| H[ll ••- H[M] 

He= H[01 H[ll H[A/] 



H[01 H[l] -.. H[M] 

is the Toeplitz diagonal-band channel matrix and: 

X [n + ij^] 



x[nl 



X [n Lf, - i\/J 
is the vector transmitted. 



(28) 
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Introducing the combined filter: 
G' = F'He - [G' [-La,] - - C [0) [Lo^]] € C^^^^^ (30) 

the following final expression is obtained: 

[n] - G'xc [n] + PV^ [n] (31 ) 

5 where Lqi = Lfi/ Lg2 = Lf2 + M and Lg = Lp + M. 

By analogy, the backward filter is defined as 
follows : 

B' - [B^ [-£bJ . . . 0 - ^ . B' [Lb,]] € C^^^^« (32) 

10 The output of the scrambling corrector may be 

written: 

[n] = G' [tjx[n ^ k] ^J^B' [k]x^[n - k] + ^F' [k] w[n-k] (33) 

k k k 

The noise and residual scrambling vector v-'-tn] (also 
called the MIMO ecjualizer error vector) is expressed as 
15 follows: 

v'[rt] - 2'[n]-G«[0]x[7i] 

- J2 G' [k] X [n » k] ^J2 b' M ^[n-k] + J2F' f^l ^ *-! 

1 

- E [* i « - - E M [n - + E ^' W ^ - ^1 (34) 

k k k 



in which: 

G' = [G' [-£oJ . . 0 . - . G' iLcl] (35) 

20 The filters and are calculated block by block 

to minimize the mean square error of the MIMO equalizer, 
vmder the constraint G^[0] = I, a problem that may be 
formulated compactly as follows in the MIMO situation: 

{F',B'} = arg luin tr E /v' [n] V [n - Arf ) 

= arg mill trfK'l 
{F3>/G<loi=i 
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The MSE minimization is executed in two successive 
steps and finally yields: 

B'-G' <37) 

P' - (EAHj*'"'HeEk)"' EaHI*'-' (38) 

where : 

0 I - ' 0 -IcC^^^^ (39) 
< s.—/^^-^,,^ f 

and where: 
with: 

(41) 

10 It should be emphasized that the noise plus 

residual scrambling at the output of the subtracter is 
correlated in space and in time. The time correlation 
has no real impact on subsequent processing but the space 
correlation has a key role in this respect. Space white 

15 noise may easily be obtained by simple Cholesky 
factorization of the correlation matrix K-*"!,. 

Proposition 2: Since the correlation matrix K^^ is 
defined as positive, Cholesky factorization is always 
applicable. Knowing this: 

Kt-LL» (42) 

20 

where L is a lower triangular matrix. Then, by applying 
F-"-' = L"-^ and ' = L-"^ B*^ as forward filter and 
backward filter (instead of and B^) , the scrambling 
plus noise correlation matrix equals the identity matrix. 
25 Proof: The proof is self-evident. 

- L-^ [{ai ^ a^^) B'B't + cr^P'P'^] (L"*)^ - L'^Kt (L^*)^ - 1 (43) 
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Validates the algorithm proposed by the Figure 3 
simulation. 

Below, the output of the equalizer is always 
considered with space whitening so that: 

z'[n] = G^[0]x[n]+v'[n] 

where Kv^ = I and G^[0] = L"^ . 

It is interesting to note that the MMSE criterion 
corresponds to maximizing the SNR tr {G^ [ 0 ] '^Kv'^G-^ [ 0 ] } . The 
maximum SNR given for the adapted filter is achieved for 



C, Exhaustive APP decoder 

All the probabilistic quantities exchanged during 
the iterations are defined below. The object of the 

15 detection portion of the MIMO equalizer is to supply 

extrinsic information on the bits of the symbols of the 
matrix A in accordance with the new communications model 
(10) where G^[0] plays the role of a flat MIMO channel 
equivalent T x t and v^[n] is the residual interference 

20 plus whitened noise. The logarithmic ratios APP over all 
the bits of the symbols a<t,j>[n] are defined as follows: 



Pr 






Pr 







(11) 
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By simple marginalization, the following is 
obtained: 



(12) 



where : 



(13) 



The logarithmic ratios APP may be expanded as 
follows : 
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,^ _ E^nO),, P ["] 1^ =t^ (a) ) P'^' [a N = aj 
Hj> N - ^ ^^^^^^^ ^ l„j |x =/x Ca) ) Pr* [a (nj = a] ^ 

Assuming that perfect space- time interleaving 
authorizes broadcasting the probabilities previously 
5 joined in the product of their marginal digit component: 

Assuming that v^[n] has been spatially whitened by 
any means (for example by Cholesky factorization), the 
Euclidean metric may be used in the evaluation of the 
10 probability: 

p (z* N |Tc) oc exp - G' (0] xf I (16) 

In the general case, the correlation matrix K"*-v of 
the noise vector must be taken into account in evaluating 
the probability: 

p (z' [n\ |x) oc exp |- (z' [n] - G' [0] x)* K^t' N - [0) x) j (17) 

As has been established in the field of 
turbodetectors, the extrinsic logarithmic probability 
rates are deduced from the equation: 

€<tj> W - n ^^^^^^ ^ j^^j 1^ J n{^^4^<i^, Pr^ [aiuj^ N ^ «{.^)] 

20 All the logarithmic extrinsic information samples 

from the MIMO detector are collected from all the blocks 
p = 1, P and rearranged, after space-time de- 
interleaving n"-'-, into a simple observation vector ^•'"^ e 
R^, on the basis whereof the output decoder supplies the 

25 logarithmic extrinsic probability ratios over all the 
code word bits: 
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Pr 



c[«i = i|g,c''7{<^''N} 



c(nl = 0|€,C''7{<'''(«l} 



(19) 



After the space- temporal de- interleaving n, the 
vector is broadcast into P matrices n^'^ of preceding 

logarithmic probability rate, one for each block of data 
A^. For each block p = 1, P, the input to n^'^ is given 
by the ecjuation: 

t' ln\ - ]n^^-tMM^lL 
Pr [a^tjf [n\ = O] 

This means that the soft decision vector x^[n] can 
be re-written as follows: 

^ w = ^ E - n (i + (-) - 0 (%^) } (21) 



p. Quasi -opt imiam MIMO detection via the sphere decoded 
modified list 

Optimum MIMO detection cannot be retained in the 

15 event of high bit rate communications scenarios because 
point enumeration of the entire subsets of constellation 
^i''°%t,j> and Q''^%t,j> # the cardinality whereof varies in 
0(2^), may rapidly become overpowering for higher order 
modulation and/or a large number of transmit antennas. 

20 Careful analysis of the probability values indicates that 
a large number of them are negligible. Consec[uently, as 
a significant contribution of this paper, we suggest 
replacing in-depth point enumeration by point enumeration 
of much smaller subsets Li^°%t,j> and L^^%t,j>/ are also 

25 called lists, which contain only non-negligible 

probabilities. The logarithmic extrinsic probability 
ratios then become: 
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In geometrical terms, they contain points of the 
trellis in a sphere of radius r centered on a carefully 
selected point (for example the point ML with no 
constraint or the point ML itself) . Modified versions of 
5 the sphere decoder have been employed with some success 
to show these lists [2] . The sphere radius r chosen 
governs the quality and complexity of the corresponding 
MIMO detector. 



10 E. PIC detection algorithm 

The sphere decoder may be considered too complex. 
One way to reduce further the complexity of the detector 
is to generate the APP log independently for each 
dimension. On iteration I, the decision variable for the 

15 component x is: 




where gk is the column k of the matrix G1[0]. The ratios 
APP log ^'^<t,j>[n], j= l,..w q are then calculated assuming 
that r^t[n] is a Gaussian variable with mean Xt and 
20 estimated variance: 




The first iteration detector may be iir^roved using 
a SIC technique in the following manner: 

• Initialization 
25 Classify the components in decreasing order as a 

function of their signal to interference ratio (SIR) , for 
example: 



SIR{xi) > SIR{x2) > ' > SIR{xt) 
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wi th: 

SIR{x,) = _ 

Recursivity 
For t = 1 to T 

5 Calculate: 

Calculate the ratio APP log ^■^<t,j>Cn] , j= 1,..., q 
assuming that r^tEn] is a Gaussian variable with mean Xt 
and estimated variance: 

1 ^ ! 2 

10 

Calculate the MMSE estimate of knowing r-^t[n]: 
xt[n] = E [xt[n] |rj[n] } 
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